Search for: All records

Understanding alpine glacier extent during past climate variability is instructive for determining the glacier response to future climate change. Villarrica volcano is a late Pleistocene stratovolcano located in Chile's Southern Volcanic Zone that was covered by the Patagonian Ice Sheet during the last glacial period, and still retains small remnant glaciers today. Moraines preserved several kilometers from the summit on different flanks of the volcano record a history of expanded glacier lengths during the Holocene. However, the precise ages of these moraines are unknown, and the Holocene glacial history of Villarrica remains poorly constrained, limiting our understanding of how glaciers in this region responded to Holocene climate change. To constrain the timing of these moraines, we analyzed cosmogenic 3He in olivine from 25 basaltic andesite moraine boulders for cosmogenic surface exposure dating. Our new chronology reveals multiple late Holocene glacier advances from different flanks of the volcano, with the glaciers culminating and abandoning their moraines during the early Neoglacial period at ∼3355 ± 190 a and ∼1735 ± 215 a, and during the last millennium spanning the Little Ice Age period at ∼720 ± 225 a, ∼370 ± 75 a, and in the last ∼200 years. Our analysis of Holocene climate proxies from south-central Chile indicates that the early Neoglacial advances and subsequent retreat likely reflect increased effective moisture delivered by intensified Southern Westerly Winds and associated shifts in their latitudinal position. In contrast, we interpret the last millennium glacier advances as primarily driven by reduced summer ablation linked to regional cooling, followed by glacier retreat due to increased temperatures. Our chronology and closely spaced moraine positions suggest that glacier retreat on Villarrica, and possibly the broader Southern Volcanic Zone, has been gradual during the late Holocene and interrupted by short-lived advances driven by varying changes in temperature and moisture. 

Abstract. Studying the retreat of the Patagonian Ice Sheet (PIS) during the last deglaciation represents an important opportunity to understand how ice sheets outside the polar regions have responded to deglacial changes in temperature and large-scale atmospheric circulation. At the northernmost extension of the PIS during the Last Glacial Maximum (LGM), the Chilean Lake District (CLD) was influenced by the southern westerly winds (SWW), which strongly modulated the hydrologic and heat budgets of the region. Despite progress in constraining the nature and timing of deglacial ice retreat across this area, considerable uncertainty in the glacial history still exists due to a lack of geologic constraints on past ice margin change. Where the glacial chronology is lacking, ice sheet models can provide important insight into our understanding of the characteristics and drivers of deglacial ice retreat. Here we apply the Ice Sheet and Sea-level System Model (ISSM) to simulate the LGM and last deglacial ice history of the PIS across the CLD at high spatial resolution (450 m). We present a transient simulation of ice margin change across the last deglaciation using climate inputs from the National Center for Atmospheric Research Community Climate System Model (CCSM3) Trace-21ka experiment. At the LGM, the simulated ice extent across the CLD agrees well with the most comprehensive reconstruction of PIS ice history (PATICE). Coincident with deglacial warming, ice retreat ensues after 19 ka, with large-scale ice retreat occurring across the CLD between 18 and 16.5 ka. By 17 ka, the northern portion of the CLD becomes ice free, and by 15 ka, ice only persists at high elevations as mountain glaciers and small ice caps. Our simulated ice history agrees well with PATICE for early deglacial ice retreat but diverges at and after 15 ka, where the geologic reconstruction suggests the persistence of an ice cap across the southern CLD until 10 ka. However, given the high uncertainty in the geologic reconstruction of the PIS across the CLD during the later deglaciation, this work emphasizes a need for improved geologic constraints on past ice margin change. While deglacial warming drove the ice retreat across this region, sensitivity tests reveal that modest variations in wintertime precipitation (∼10 %) can modulate the pacing of ice retreat by up to 2 ka, which has implications when comparing simulated outputs of ice margin change to geologic reconstructions. While we find that TraCE-21ka simulates large-scale changes in the SWW across the CLD that are consistent with regional paleoclimate reconstructions, the magnitude of the simulated precipitation changes is smaller than what is found in proxy records. From our sensitivity analysis, we can deduce that larger anomalies in precipitation, as found in paleoclimate proxies, may have had a large impact on modulating the magnitude and timing of deglacial ice retreat. This fact highlights an additional need for better constraints on the deglacial change in strength, position, and extent of the SWW as it relates to understanding the drivers of deglacial PIS behavior. 

Tropical glaciers have retreated over recent decades, but whether the magnitude of this retreat exceeds the bounds of Holocene fluctuations is unclear. We measured cosmogenic beryllium-10 and carbon-14 concentrations in recently exposed bedrock at the margin of four glaciers spanning the tropical Andes to reconstruct their past extents relative to today. Nuclide concentrations are near zero in almost all samples, suggesting that these locations were never exposed during the Holocene. Our data imply that many glaciers in the tropics are probably now smaller than they have been in at least 11,700 years, making the tropics the first large region where this milestone has been documented. 

The last glacial period was punctuated by cold intervals in the North Atlantic region that culminated in extensive iceberg discharge events. These cold intervals, known as Heinrich Stadials, are associated with abrupt climate shifts worldwide. Here, we present CO₂measurements from the West Antarctic Ice Sheet Divide ice core across Heinrich Stadials 2 to 5 at decadal-scale resolution. Our results reveal multi-decadal-scale jumps in atmospheric CO₂concentrations within each Heinrich Stadial. The largest magnitude of change (14.0 ± 0.8 ppm within 55 ± 10 y) occurred during Heinrich Stadial 4. Abrupt rises in atmospheric CO₂are concurrent with jumps in atmospheric CH₄and abrupt changes in the water isotopologs in multiple Antarctic ice cores, the latter of which suggest rapid warming of both Antarctica and Southern Ocean vapor source regions. The synchroneity of these rapid shifts points to wind-driven upwelling of relatively warm, carbon-rich waters in the Southern Ocean, likely linked to a poleward intensification of the Southern Hemisphere westerly winds. Using an isotope-enabled atmospheric circulation model, we show that observed changes in Antarctic water isotopologs can be explained by abrupt and widespread Southern Ocean warming. Our work presents evidence for a multi-decadal- to century-scale response of the Southern Ocean to changes in atmospheric circulation, demonstrating the potential for dynamic changes in Southern Ocean biogeochemistry and circulation on human timescales. Furthermore, it suggests that anthropogenic CO₂uptake in the Southern Ocean may weaken with poleward strengthening westerlies today and into the future. 

Abstract. Determining the timing and extent of Quaternary glaciations around the globe is critical to understanding the drivers behind climate change and glacier fluctuations. Evidence from the southern mid-latitudes indicates that local glacial maxima preceded the global Last Glacial Maximum (LGM), implying that feedbacks in the climate system or ice dynamics played a role beyond the underlying orbital forcings. To shed light on these processes, we investigated the glacial landforms shaped and deposited by the Lago Argentino glacier (50° S), an outlet lobe of the former Patagonian Ice Sheet, in southern Argentina. We mapped geomorphological features on the landscape and dated moraine boulders and outwash sediments using 10Be cosmogenic nuclides and feldspar infrared stimulated luminescence (IRSL) to constrain the chronology of glacial advance and retreat. We report that the Lago Argentino glacier lobe reached more extensive limits prior to the global LGM, advancing during the middle to late Pleistocene between 243–132 ka and during Marine Isotope Stage 3 (MIS 3), culminating at 44.5 ± 8.0 and at 36.6 ± 1.0 ka. Our results indicate that the most extensive advance of the last glacial cycle occurred during MIS 3, and we hypothesize that this was a result of longer and colder winters, as well as increased precipitation delivered by a latitudinal migration of the Southern Westerly Winds belt, highlighting the role of local and regional climate feedbacks in modulating ice mass changes in the southern mid-latitudes. 

Abstract Subannual climate reconstructions of the Holocene are rare despite the ability of such records to provide a better understanding of the underlying factors that drive subannual climate variability. We used specialized confocal laser fluorescent microscope imaging and automated secondary ion mass spectrometry microanalysis to resolve a seasonal oxygen isotope (δ18O) record of a late Holocene–aged (2.7–2.1 ka) speleothem from mid-continental North America. We did this by measuring intra-band δ18O variability (Δ18O) within 117 annual bands over a 600 yr span of the late Holocene. We interpret a change in Δ18O values after 2.4 ± 0.1 ka to reflect an increase in the amount of winter precipitation. Our study produced direct measurements of past seasonality, offers new insights into shifting seasonal precipitation patterns that occurred during the late Holocene in central North America, and adds a new tool for understanding the complex precipitation and temperature histories of this region. 

Abstract. There is unambiguous evidence that glaciers have retreated from their 19th century positions, but it is less clear how far glaciers have retreated relative to their long-term Holocene fluctuations. Glaciers in western North America are thought to have advanced from minimum positions in the Early Holocene to maximum positions in the Late Holocene. We assess when four North American glaciers, located between 38–60∘ N, were larger or smaller than their modern (2018–2020 CE) positions during the Holocene. We measured 26 paired cosmogenic in situ 14C and 10Be concentrations in recently exposed proglacial bedrock and applied a Monte Carlo forward model to reconstruct plausible bedrock exposure–burial histories. We find that these glaciers advanced past their modern positions thousands of years apart in the Holocene: a glacier in the Juneau Icefield (BC, Canada) at ∼2 ka, Kokanee Glacier (BC, Canada) at ∼6 ka, and Mammoth Glacier (WY, USA) at ∼1 ka; the fourth glacier, Conness Glacier (CA, USA), was likely larger than its modern position for the duration of the Holocene until present. The disparate Holocene exposure–burial histories are at odds with expectations of similar glacier histories given the presumed shared climate forcings of decreasing Northern Hemisphere summer insolation through the Holocene followed by global greenhouse gas forcing in the industrial era. We hypothesize that the range in histories is the result of unequal amounts of modern retreat relative to each glacier's Holocene maximum position, rather than asynchronous Holocene advance histories. We explore the influence of glacier hypsometry and response time on glacier retreat in the industrial era as a potential cause of the non-uniform burial durations. We also report mean abrasion rates at three of the four glaciers: Juneau Icefield Glacier (0.3±0.3 mm yr−1), Kokanee Glacier (0.04±0.03 mm yr−1), and Mammoth Glacier (0.2±0.2 mm yr−1).